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The  paper  continues  a  series  of  publications  devoted  to  modem  advances  in  aerodynamic  decelerator  system 
technology  started  recently  ( Journal  of  Aircraft ,  Vol.  38,  No.  5,  2001)  and  addresses  the  development  of  a  six- 
degree-of-freedom  model  of  a  guided  circular  parachute.  The  paper  reviews  existing  circular  parachute  models 
and  discusses  several  modeling  issues  unresolved  within  the  frame  of  existing  approaches  or  completely  ignored  so 
far.  These  issues  include  using  data  obtained  in  the  aerodynamic  experiments  and  computational-fluid-dynamics 
modeling  for  both  undistorted  (uncontrolled)  and  distorted  (controlled)  canopy  shapes,  introducing  and  computing 
control  derivatives,  and  providing  comparison  with  the  real  flight  data.  The  paper  provides  step-by-step  develop¬ 
ment  of  the  mathematical  model  of  circular  parachute  that  includes  the  basic  equations  of  motion,  analysis  and 
computation  of  the  aerodynamic  forces  and  moments,  and  investigation  with  modeling  of  special  modes  observed 
in  flight.  It  then  introduces  a  new  application  of  a  two-step  aerodynamic  parameters  identification  algorithm  that 
is  based  on  comparison  with  two  types  of  the  air-drop  data  (uncontrolled  set  and  controlled  one).  The  paper  ends 
with  summary  of  the  obtained  results  and  proposes  a  vital  direction  for  the  further  elaboration  of  the  developed 
model. 
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Introduction 


PARACHUTES  have  been  the  simplest  and  cheapest  devices 
used  for  the  delivery  of  materials,  people,  and  vehicles  ever 
since  their  first  recorded  use  by  Jacques  Garnering  who  jumped 
from  a  balloon  over  Paris  in  1797.  However,  this  very  simplicity 
makes  their  aerodynamics  very  difficult  to  model.  Specifically,  the 
parachute  not  only  deflects  surroundingair,  but  also  adopts  its  shape, 
which  is  dictated  by  the  airflow  generated  by  the  canopy.  Further¬ 
more,  lack  of  streamliningduring  controlactivationmakes  turbulent 
rather  than  laminar  flow  dominate  the  parachute’s  aerodynamics. 

Significant  research  on  flat  circular  parachute  modeling  has  been 
done  over  the  past  60  years  by  researchers  in  the  U.S.,  Europe, 
and  Russia.  However,  no  complete  model  of  a  controlled  circular 
parachute  has  been  developed  so  far  (whereas  a  number  of  efficient 
techniques  addressing  the  control  of  maneuverable  parachutes  like 
parafoil  exists).  Moreover,  existing  high-degree-of-freecbm  models 
of  circularnoncontrolledparachuteslack  verified  nonlinearaerody- 
namicsand  use  mostly  empirical  values  for  the  apparentmass  terms 
(see  Refs.  1-22). 

The  main  contribution  of  the  work  reported  here  is  the  devel¬ 
opment  of  a  controlled  model  of  a  six-degree-of-freedcm  (DoF) 
circular  parachute.  This  was  done  by  applying  a  nonlinear  system 
identification  algorithm  to  refine  the  computational-fluid-dynamics 
(CFD)  values  of  the  aerodynamic  coefficients.  The  critical  idea  was 
to  first  identify  the  uncontrolled  model  aerodynamics  and  use  the 
resulting  estimates  to  identify  the  controlled  model  aerodynamics. 
Another  important  step  was  to  subtract  the  influence  of  the  wind 
when  determining  the  controlled  model.  This  approach  was  made 
possible  by  the  richness  of  the  flight  data  set,23'24  that  is,  data  for 
both  uncontrolled  and  controlled  drops  were  available. 

A  detailed  analysis  of  the  circular  parachute  modeling  problem 
discussed  in  the  literature  has  identified  two  shortcomings  inherent 
in  existing  methods:  1)  estimation  of  apparentmass  terms  is  done 
empirically  for  axisymmetric  shapes;  2)  the  only  aerodynamics  con¬ 
sidered  are  those  of  a  fully  deployed  and  symmetric  canopy. 

Based  on  Henn’s1  work  introducing  apparent  mass  coefficients 
(1944),  several  authors  developed  mathematical  models  of  a  dif¬ 
ferent  degree  of  complexity  to  investigate  the  dynamic  behavior  of 
parachutes  (e.g.,  see  Ref.  2).  Different  empirical  values  for  two  ap¬ 
parent  masses  were  used  in  these  studies.  In  1962  Lester3  rederived 
basic  equations  of  motion  showing  that  Henn’s  equations  were  er¬ 
roneous.  At  the  same  time  Ibrahim4-5  has  conducted  the  first  ex¬ 
haustive  analysis  of  the  apparent  mass  and  moment  of  inertia  terms 
of  cup-shapedbodies  in  unsteady  incompressible  flow.  He  was  also 
the  first  one  who  showed  that  resulting  parachute  dynamic  perfor¬ 
mance  is  sensitive  to  the  values  of  these  terms.  This  study  initiated 
the  scientific  debate  on  the  extent  of  influence  of  the  apparentmass 
terms  on  the  stability  of  a  descending  system.  Many  authors13-17 
concluded  that  neglecting  nonlinear  terms  and  the  stochastic  nature 
of  parachute  dynamics  would  lead  to  completely  inaccurate  model¬ 
ing  and  also  stressed  the  need  for  experimental  determination  of  the 
apparent  mass  terms.  An  attempt  to  use  theoretical  apparent  mass 
and  moments-of-inertiacoefficients  based  on  those  related  to  the  air 
displacedby  ellipsoids  of  revolution  moving  in  a  potential  fluid  was 
made  also.12  One  of  the  latest  experimental  studies14' 15  suggests  that 
apparentmasses  and  moments  of  inertia  (hereafterreferred  to  jointly 
as  the  apparentmass  tensor  terms)  are  the  functionsof  a  spatial  angle 
of  attack  and  depending  on  acceleration  might  significantly  exceed 
the  theoretical  ones. 


With  regard  to  aerodynamics  of  the  fully  deployed  canopies  with 
a  symmetric  shape,  various  authors  have  shown  that  the  nonlin¬ 
ear  nature  of  basic  aerodynamic  terms  is  a  function  of  the  angle 
of  attack.  For  instance,  in  Ref.  16  authors  discuss  parameter  iden¬ 
tification  using  flight-test  data,  where  they  point  out  that  the  de¬ 
pendence  of  aerodynamic  terms  on  the  angle  of  attack  is  highly 
nonlinear  unlike  those  of  an  aircraft.  Many  of  papers12,14,15,17  in¬ 
clude  sufficient  data  for  fully  deployed  canopies  and  a  variety  of 
flight  conditions.  (Similar  to  the  aircraft  aerodynamics  the  major¬ 
ity  of  available  data  has  been  obtained  on  the  base  of  wind-tunnel 
experiments. )  However,  the  results  reported  employ  predetermined 
symmetric  canopy  shapes.  The  aerodynamics  of  a  distorted  canopy 
has  not  been  considered. 

In  summary,  over  the  past  40  years  a  consensus  on  the  lack  of 
accurate  dynamic  modeling  of  apparent-mass  effects  and  nonlinear 
aerodynamics  of  distorted  canopies  has  emerged. 

The  complexity  of  parachute  motion  has  also  been  confirmed  by 
the  extensive  amount  of  real  air-drop  experiments.21-24  First,  it  has 
been  verified  that  strongly  non-linear  nature  of  the  parachute  aero¬ 
dynamics  is  determined  by  three  major  factors:  flight  conditions, 
canopy  shape,  and  geometry  of  “parachute-payload”system  [there¬ 
after  referred  to  as  an  aerial  delivery  system  (ADS)] .  Thus ,  assuming 
that  the  material  geometry  does  not  change  during  flight  leads  to  the 
conclusion  that  aerodynamic  forces  and  moments  depend  on  canopy 
shape,  spatial  angle  of  attack,  and  dynamic  pressure. 

Second,  the  distortionof  the  axisymmetric  canopy  shape  rules  out 
employing  the  well-established  analytical  results  used  to  determine 
the  apparentmass  tensor  for  axisymmetric  canopy.  This  problemhas 
not  been  resolved.  Therefore,  based  on  preliminary  results  provided 
by  CFD  techniques25  and  our  own  analysis,  several  assumptions 
have  been  made  to  model  canopy  distortioncausedby  control  inputs. 
In  particular,  it  was  assumed  that  the  number  of  canopy  shapes 
caused  by  control  inputs  are  finite  and  known. 

At  the  same  time  there  is  an  essential  yaw  rotation  caused  by 
canopy  shape  distortion  when  even  a  small  difference  between  the 
length  of  two  adjacent  control  inputs  is  presented.  Therefore,  the 
control  input  performance  strongly  affects  the  whole  system  dy¬ 
namics.  Hence,  the  controlled  parachute  model  should  properly 
consider  it. 

Finally,  the  descentrate  dependsnot  only  on  air  density  but  also  on 
control  inputs  through  the  canopy  shape  distortion  that  was  also  ob¬ 
served  during  the  real  drops.  Therefore,  the  model  of  the  controlled 
flat  circular  parachute  should  include  six  degrees  of  freedom  with 
an  appropriate  description  of  apparent  mass  tensor,  aerodynamics, 
and  control  system  dynamics. 

This  paper  is  organizedas  follows.  The  next  section  addresses  the 
development  of  a  six-DoF  parachute  model.  It  contains  the  detailed 
discussionof  the  mathematical  issues  involved,  beginning  with  anal¬ 
ysis  of  the  basic  equations  of  motion,  followed  by  computation  of 
aerodynamic  forces  and  moments,  and  ending  with  the  modeling  of 
some  special  modes  observed  in  the  flight  test.  The  following  sec¬ 
tion  discusses  verification  of  the  parachute  model  using  a  two-step 
parameter  identification  technique.  The  paper  ends  with  the  conclu¬ 
sions.  The  Appendix  contains  mass-geometry  datum  for  a  generic 
ADS  consisting  of  a  G- 12  parachute  and  an  A-22  container,  which 
was  used  for  modeling. 

Model  Development 

Definitions  and  Assumptions 

In  this  section  we  discuss  the  development  of  a  circularparachute 
model.  The  specific  parameters  and  geometry  of  descending  sys¬ 
tem  used  were  those  of  a  G-12  parachute.  The  G-12  is  a  150-m2 
nylon  cargo  parachute  with  64  suspension  lines  (SL)26  (Fig.  1 ).  The 
A-22  delivery  container  was  selected  as  a  prototype  cargo  box.  This 
1.82-m3  almost  cubic  container  has  a  payload  capacity  of  nearly  a 
ton  and  is  commonly  used  with  the  G-12  by  the  U.S.  Army. 

The  following  controlled  ADS  architecture  is  considered.  All 
lines  are  assembled  into  eight  link  assembles  (see  Fig.  1).  Each 
pair  of  assembles  is  attached  to  one  of  four  risers.  At  the  other  end 
the  risers  are  coupled  to  the  payload  at  four  dispersed  points  (see 
Fig.  2).  By  design  these  risers  allow  controlling  the  parachute  by 
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Fig.  3  System  of  axes  and  definition  of  aerodynamic  terms. 


lengthening  one  or  two  adjacent  actuators  and  hence  disturbing  the 
symmetric  shape  of  the  canopy.  The  pneumatic  muscle  actuators 
(PMA)  developed  by  Vertigo,  Inc.,27  are  modeled  and  used  in  this 
study  as  a  prototype  of  the  control  risers. 

Two  Cartesian  coordinate  systems  have  been  chosen  to  describe 
ADS  motion  (see  Fig.  3). 

Linear  position  of  the  ADS  is  computed  with  respect  to  a  local 
tangent  plane  {/<}.  Its  positive  y  directionis  aligned  with  local  North, 
the  positive*  directi onpoints  East,  and  the  positive;  directi onpoints 
up. 

All  other  computations  are  performed  in  the  body-fixed  coordi¬ 
nate  frame  {£>} .  Its  origin  is  attached  to  the  center  of  the  open-end 
plane  of  the  canopy.  The  x  and  y  body  axes  lie  in  the  plane  paral¬ 
lel  to  the  canopy’s  base,  and  z  is  aligned  with  the  imaginary  axis 
extending  toward  the  centroid  of  the  payload. 

It  is  worth  mentioning  thatin  differentaerohydrodynamicstudies 
the  origin  of  {b}  is  sometimes  placed  at  the  canopy’s  centroid.  The 
undisturbed  canopy  is  always  assumed  to  be  a  planetary  ellipsoid, 
but  it  can  have  a  differentratio  of  minor  and  major  axes.  That  makes 
the  location  of  the  origin  of  {i>}  to  be  conditional  from  the  concrete 
parachute  design.  Moreoverthe  z  coordinateof  the  canopy  centroid 
can  be  determined  only  approximately  because  in  the  general  case 
there  is  no  analytical  formula  for  ellipsoidal  shells.  Despite  the  fact 
that  originally  the  six-DoF  model  presented  in  this  paper  was  also 
developed  with  this  setup  of  the  frame  {£>},  it  was  later  on  redone 
with  the  origin  at  the  center  of  the  open-end  plane  of  the  canopy  to 
make  equations  of  motion  more  universal. 

For  simplicity  we  will  skip  subscript  b  in  the  further  analysis, 
assuming  that  all  variables  and  aerodynamic  coefficients  when  ap¬ 
plicable  are  defined  in  {£>} . 

The  following  assumptions  mostly  adopted  from  the  Tory  and 
Ayres  paper12  were  used  to  develop  the  model: 

1 )  Because  of  the  predetermined  architecture,  the  parachute  and 
payload  are  consideredto  be  rigidly  connected  to  each  other  (Fig.  2). 

2)  During  the  air  drop,  these  two  rigidly  connected  parts  are  as¬ 
sumed  to  experience  only  gravity  and  aerodynamic  forces. 

3)  The  canopy  experiences  all  aerodynamic  forces  and  moments 
about  its  center  of  pressure. 


4)  The  aerodynamic  forces  generated  by  the  payload  are  negligi¬ 
ble. 

5)  Undistorted  canopy  has  an  axial  symmetry  about  the  z  axis. 

Distortion  of  the  canopy’s  shape  (causedby  the  lengtheningof  one 

or  two  adjacent  risers)  introduces  asymmetric  forces  and  moments 
allowing  steering  of  the  ADS  in  a  certain  direction  in  a  horizontal 
plane.  Moreover,  it  obviously  makes  ADS  nonsymmetrical.  (One 
plane  of  symmetry  still  remains,  but  its  location  is  not  constant  with 
respectto  the  body  frame.)  In  this  study  however  we  will  additionally 
assume  that  the  effect  of  risers  lengthening  and  canopy  distortion 
from  the  standpoint  of  changing  tensors  of  inertia  and  apparent 
masses  is  negligibly  small.  Therefore  even  for  a  controlled  circular 
parachute  these  two  tensors  will  be  assumed  to  have  the  same  form 
as  for  uncontrolled  (symmetrical)  parachute. 

Aerial  Delivery  System  Geometry 

Figure  4  introduces  a  set  of  parameters  defining  the  parachute’s 
geometry.  Numbers  1 ...  4  denote  the  ADS’s  components:  canopy, 
rigging  (suspension)lines,  actuators  (PMAs),  and  payload.  The  po¬ 
sition  of  their  centroids  z,  are  determined  with  respect  to  the  body 
frame  origin  (point  O  on  Fig.  4). 

On  this  figure  Rp  =  jR0  denotes  the  radius  of  inflated  parachute, 
and  Opi  denotes  the  dimension  of  a  container  (without  loss  of  gen¬ 
erality  all  three  dimensions  of  a  cargo  container  were  consideredto 
be  the  same). 

As  just  mentioned,  the  shape  of  undisturbed  canopy  is  a  half  of 
a  planetary  ellipsoid  (hemispheroid),  meaning  that  it  is  circular  in 
plan  when  viewed  along  the  z  direction  and  elliptical  when  viewed 
from  the  side.  The  ratio  of  the  minor  to  major  axes  (canopy  shape 
ratio)  will  be  denoted  as  e.  Tory  and  Ayres12  reported  that  this  ratio 
for  a  flat  circular  canopy  is  typically  equal  to  0.5.  Other  researchers 
(e.g.,  see  Refs.  14-17)  assumed  canopy  to  be  a  hemisphere  (s  =  1). 
For  the  ADS  at  hand,  e  is  equal  to  0.82. 

Basic  Six-Degree-of-Freedom  Equations 

Following  basic  analytical  mechanics  principles,  the  total  energy 
T  of  a  whole  descending  system  relative  to  body  frame  can  be 
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given  by 

y  _  yADS  _j_  yair  /j\ 

The  expressions  for  T ,  TADS,  and  Tml  can  be  found  in  Ref.  16. 

For  an  ideal  fluid  the  kinetic  energy  can  be  determinedin  terms  of 
velocity  potential.  However,  it  is  common  practice16  to  assume  that 
the  kinetic  energy  of  real  fluid  can  be  defined  similarly.  Using  the 
Lagrange  approach,  the  basic  equations  of  motion  for  the  parachute- 
air  system  can  be  obtained.  The  most  general  form  of  these  equations 
in  (b)  is 


d 

dt 


+  ST  x 


d 

dr 


+  ST  x 


+  V  x 


(2) 

(3) 


where  after  substitutingthe  expressionfor  kinetic  energy  of  the  body 
and  the  air  into  Eqs.  (2)  and  (3)  and  applying  appropriate  apparent 
mass  tensor  for  a  symmetric  body,  and  some  algebra,  the  final  form 
of  equations  of  motion  will  appear  to  be  as  follows: 


rigging  lines  m2,  actuators  m 3,  and  payload  m4 
4 

K  =  y  nijZj  =  mzG 
i  =  1 

where  zG  is  the  static  c.g.  of  the  overall  ADS  with  respect  to  the  point 
O  (see  Fig.  4).  The  other  notations  denote  the  diagonal  components 
of  ADS’s  inertia  tensor  Ixx,  Iyy,  j /„  and  apparentmass  tensor's  com¬ 
ponents  amn.  The  only  asymmetry  left  in  these  equations  is  caused 
by  a  possible  asymmetry  of  payload.  (Ixx  in  general  might  not  be 
equal  to  Iyy.)  (Computation  of  both  tensors  is  considered  in  details 
in  the  following  subsections.) 

Equations  (4)  and  (5)  are  very  similar  to  those  used  to  model  a 
rigid-body  motion  [F  =  m  V  +  m (ST  x  V )  and  M  =  /ST  +  ST  x  /ST], 
and  in  vector  form  they  can  be  rewritten  as  follows: 

F  =  Mm  V  +  A cMm ,  M  =  /ST  +  HI  +  Mcr  (6) 

In  the  latter  expression 


Mm  =  diagjm  +  an,m  +  an ,  m  +  a33) 
— vr  wq  q  +  rp 

Ac  =  ur  —wp  qr  —  p 

_vp  —  uq  0  —  (p2  +  q2)  J 

Mm  =  [m  +  an,  m  +  a33,  K  +  a15]r 
I  =  diag{/v  v  -{-  044,  Iyy  -f-  a44,  Iz:  -F  a66} 


H,  = 


wp  —  ur  —  v  —qr  vw 
wq  —  vr  +  u  pr  —uw 
0  0  0 


I  —  [ K  +  an,  Iyy  +  a44  —  I:z  —  a66,  a33  —  an  J 

Mcr  =  L0,  0,  (Iyy  -  Ixx)pq]T 


Being  resolved  with  respect  to  V  and  ST,  Eqs.  (6)  yield 

V  =  M~'(F- AcMm),  ST  =  I~l(M  —  HCI  —  Mcr)  (7) 

The  attitude  of  the  ADS  is  determined  by  the  Euler  angles  (p,  1}, 
and  i/r .  The  relation  between  vector  ST  and  Euler  angles  is  found  in 
the  usual  way28,29: 


<i> 

"1 

sin  ip  tg  8 

cos  <j>  tg  6 

e 

- 

0 

COS  (f> 

sin  (j> 

-i- 

_o 

sin0(l/cos@) 

cos  0(1 /cost?) 

The  local  tangent  plane  coordinatesP,,  =  [x„,  y„ ,  z„]T  of  the  ori¬ 
gin  of  {/>}  can  be  obtained  employing  correspondingrotation  matrix 
“P: 

P  u  =  "„RV  (9) 


F  = 


(in  +  an)(fi  —  vr)  +  (m  +  a33)wg  +  (K  +  ai5)(^  +  rp) 
(m  +  an)(u  +  ur)  —  (m  +  a33)wp  —  (K  +  ai5)(/)  —  qr) 
( m  +a33)u)  —  (m  +  an)(uq  —  vp)  —  (K  +ai5)(p2  +  q2) 


M  = 


(Ixx  +  a44)p  -  (K  +  al5)(v  -  wp  +  ur)  -  (Iyy  +  a44  -  /„  -  a66)qr  +  (a33  -  an)vw 
(Iyy  +  u44)q  +  (K  +  a15)(u  +  wq  -  vr)  +  (Iyy  +  a44  -  /„  -  u66)pr  -  (a33  -  an)#w 
(Fz  +  a66)f  +  (Iyy  -  Ixx)pq 


(4) 


(5) 


(An  example  of  step-by-step  derivation  of  these  equations  can  be 
found  for  example  in  Ref.  16.) 

In  Eqs.  (4)  and  (5), 


4 


i=  1 


denotes  the  mass  of  the  system  that  includes  the  mass  of  canopy  m  { , 


Computation  of  Moments  of  Inertia 

The  static  mass  center  and  moments  of  inertia  are  determined 
based  on  the  weight  and  dimensions  of  each  component. 

First,  the  z  coordinate  of  each  component  centroid  with  respect 
to  one  of  its  surfaces  z,  and  individual  central  moments  of  inertia 
for  each  component  in  the  axes  of  correspondent  centroid  /  '  .  were 
derived.  (Magnitudes  of  z,  then  can  be  computed  with  account  of  z; .) 


486 


DOBROKHODOV,  YAKIMENKO,  AND  JUNGE 


Table  1  Relevant  formulas  of  moments  of  inertia  for  parachute  components 


Component 

Canopy  (1) 

Suspension  lines  (2)  and  PMAs  (3) 

Payload  (4) 

Z;  =  |B?| 

RP  0.83 

2 

L  cos  y 

2 

Clp\ 

T 

=  4  b 

0.248m  iR2pe°-S2e 

0.246 mlR2pe°(lfx  =  lfy) 

m* 

T 

L2(l+  cos2  y) 

12 

,n4api 

6 

4  c 

^■mi«2(l +  0.143  la  s') 

(  L 2  sin2  y  9  \ 

H  12  +R*) 

m4a2 

6 

z=c 

A 


z=c 


z=c 


Table  1  contains  formulas  for  each  component.  Figures  5a-5c 
show  the  components.  Although  formulas  for  a  hemispherical  shell 
and  a  solid  cube  (payload)  are  well  known  (for  example,  see  Ref .  30), 
the  formulas  for  a  hemispheroidal  shell  and  for  a  frustum  right  cir¬ 
cular  cone  shell  formed  by  SLs  or  PMAs  are  original.  Approximate 
formulas  for  a  hemispheroidal  shell  were  specifically  derived  to 
match  real  values  in  the  range  of  e  e  [0.5;  1],  (Exact  cumbersome 
formulas  for  a  hemispheroidal  shell  can  be  found  for  instance  in 
Ref.  31.) 

In  Table  1  Rt  is  a  radius  of  the  shell  measured  at  the  z  co¬ 
ordinate  of  its  centroid.  (R„  =  Rp  —  /SL  sin  y/2  for  SLs  and  R,  = 
cZpi /  \J2  T  /pma  sin  y/2  for  PMAs.)  The  cone  half-angle  y  for  the 
considered  ADS  configuration  can  be  computed  from  the  geometric 
relation  sin  y  =  (Rp  -  apl/ s/2)(lSL  +  /pma)-1. 

The  individual  inertia  components  were  then  transferred  to  the 
origin  of  {/>}  using  the  parallel  axis  theorem 


Finally,  the  moments  of  inertia  for  the  whole  ADS  were  computed 
as  a  sum  of  inertias  of  corresponding  ADS  components 

4 

/..  •  V  /.. 

11  /  i  jj 

i  =  1 

The  numerical  values  for  the  moments  of  inertia  of  each  compo¬ 
nent  of  a  generic  ADS  are  given  in  Table  A4  in  the  Appendix.  Here 
it  can  be  stated  that  although  the  major  contributions  into  Ixx  and 
Iyy  moments  of  inertia  are  obviously  caused  by  the  payload  (more 
than  98%),  moment  Izz  is  formed  basically  by  the  parachute  itself 
(canopy,  SLs,  and  PMAs)  with  a  minor  influence  of  payload  (less 
than  15%).  That  shows  why  it  is  reasonable  to  neglect  the  effect  of 
ADS  asymmetry  while  lengthening  risers.  Another  feature  is  that 
the  symmetry  of  the  cargo  box  with  respect  to  the  axis  z  simplifies 
Eqs.  (7)  zeroing  vector  Mcr. 

Apparent  Mass  Terms 

Similar  to  the  rigid-body  mass  tensor,  the  apparent  (virtual)  mass 
tensor  A  has  6-6  =  36  elements,  and  in  a  real  fluid  these  can  all 
be  unique.  For  an  ideal  fluid,  however,  A  is  a  symmetrical  tensor, 
leaving  a  maximum  of  21  distinct  terms.  In  the  case  of  a  body 
with  two  planes  of  symmetry  and  coordinate  frame  origin  located 
somewhere  on  the  axis  of  symmetry,  tensor  A  can  be  furtherreduced 


to  the  following  form: 

"an 

0 

0 

0 

£*15 

0  ■ 

0 

a22 

0 

£*24 

0 

0 

0 

0 

a33 

0 

0 

0 

A  = 

0 

0 

0 

(ID 

a24 

0 

£*44 

a  15 

0 

0 

0 

£*55 

0 

.  0 

0 

0 

0 

0 

t^66— 

Here  the  first  three  diagonal  elements  represent  apparent  masses 
of  the  air  virtually  stagnant  within  and  around  (below)  the  canopy, 
the  next  three  are  correspondent  apparent  moments  of  inertia,  and 
the  off-diagonal  air  mass/inertia  elements  contributeto  the  coupling 
motion. 

Because  of  axial  symmetry  of  the  circular  canopy,  £*22  =  £*iit 
of 55  =  a44,  anda24  =  —  aq5.  That  leaves  only  three  distinct  elements, 
which  are  oq  i ,  c*33 ,  a44 ,  oq5 ,  and  c*66 .  [That  is  why  dynamic  equations 
have  their  appearance,  Eqs.  (4)  and  (5)J. 

In  the  earlier  studies  to  represent  a  flow  around  a  fully  deployed 
canopy,  the  latter  was  represented  as  a  spheroid.  In  this  case  the 
reference  air  mass  and  moments  of  inertia  correspond  to  those  of 
the  air  displaced  by  the  body 

ma  =  m*  =  ^npRps  (12) 

C  =  \maR\{X  +  s2),  Iff  =  f  maR\  (13) 
(see  conventions  of  Table  1). 

Today  it  is  a  usual  practice  to  refer  to  the  air  trapped  within  a 
hemisperoid.  In  this  case  air  mass  makes  the  half  of  spheroid 

ma  =  0.5m*  (14) 


Ih  =  IJJ  +  m>zi 


and  formulas  for  the  moments  of  inertia  are  the  same  as  Eqs.  (13).  If 
needed  to  be  computed  with  respectto  centroid  axes,  the  z  coordinate 


i  =  1 . 4 
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for  the  hemispheroid's  centroid  equal  to  zP  —  —  \eRp  needs  to  be 
taken  into  consideration. 

Generally  speaking,  apparent  mass  terms  depend  on  the  canopy’s 
configuration,  porosity,  acceleration,  and  spatial  angle  of  attack.  For 
instance,  Ibrahim4'5  showed  that  apparent  masses  can  drop  their  val¬ 
ues  more  than  20  times  with  the  porosity  increased  from  0  to  40%. 
Apparent  moments  of  inertia  also  decrease  their  values  by  the  fac¬ 
tor  of  2.75.  In  their  experimental  study  Yavuz  and  Cockrell14  and 
Yavuz15  demonstrated  strong  dependence  of  the  angle  of  attack  on 
especially  a33.  (It  decreases  4.5  times  with  increase  of  the  angle  of 
attack  from  0  to  40  deg.)  They  also  revealed  a  strong  dependence 
of  apparent  mass  terms  from  appropriate  accelerations  showing  that 
they  can  change  as  much  as  by  a  factor  of  five  while  experiencing 
steady  acceleration.  With  regard  to  the  latter  study,  however  it  is 
unclear  how  the  dependence  from  acceleration  can  be  taken  into 
account.  First,  the  parachute  does  not  experience  a  constant  accel¬ 
eration  during  descent.  According  to  the  flight-test  data  available, 
accelerations  during  descent  of  a  fully  deployed  parachute  oscillate 
around  zero.  But  even  the  value  at  zero  acceleration  cannot  be  im¬ 
plemented  because  all  dependences  have  a  first-order  discontinuity 
at  this  point  so  that  the  values  differ  as  much  as  three  times  when 
approaching  to  the  steady-state  descent  (zero  acceleration)  from 
negative  to  positive  acceleration,  not  mentioning  that  the  data  were 
obtainedfor  the  small  rigid  models  of  a  hemisphere  canopy.  Second, 
the  physical  values  of  apparent  masses  obtained  at  zero  acceleration 
when  being  scaled  to  the  whole  ADS  make  no  sense  because  they 
exceed  those  of  ADS  itself  (see  the  following). 

Therefore  in  the  present  study  the  authors  followed  all  other  major 
studies  and  considered  all  apparent  mass  terms  to  depend  explicitly 
on  the  air  density  only.  All  other  possible  effects  were  represented 
by  constant  coefficients: 

®n  =  knina,  <233  =  k33ma,  a44  =  k44lxxa 

0^66  —  0^15  —  k\3fflaZp  (15) 

As  opposed  to  other  studies17  where  expression  lB  — 
,/[(/SL  +  W)2  ~  R2]  (see  Fig.  4)  was  used  as  a  reference  length, 
zP  representing  the  distance  from  frame  {£>}  origin  to  the  point  of 
application  of  the  translationalapparentmass  component  (canopy’s 
center  of  pressure)  was  used  because  of  the  sense  of  Eqs.  (4)  and 

/g\  4,5,9-11,14 

Following  Ref.  17,  coefficients  in  Eqs.  (15)  were  taken  to  be  as 
follows:  kn  =0.5,  fc33  =  1.0,  k44  =0.24,  and  kl5  =0.75.  [Actually 
the  values  given  are  slightly  different  from  the  values  in  Ref.  17 
because  they  were  recomputed  to  match  notations  (13)  and  (14).] 

For  the  ADS  at  hand  with  ma  =  472  kg  (p  =  1  kgm-3),  the  follow¬ 
ing  expressions  are  valid:  an  =  236  kg  (~22%  of  m),  a33  =  472  kg 
(~44%  of  m),  a44  =  1600  kgm2  (less  than  0.5%  of  Ixx),  and 
aq5  =707  kgm  (~3%  of  K).  By  looking  at  these  numbers,  it  be¬ 
comes  clear  that  an  and  a33  are  evidently  of  most  importance.  The 
difference  between  them  affects  the  dynamic  behaviorof  the  system 
[see  Eqs.  (5)J. 

To  give  an  impression  what  coefficients  were  used  by  the  other 
researches.  Table  2  contains  some  data  on  the  major  research  in  the 
area  for  the  past  60  years. 

Henn1  was  the  first  to  perform  a  stability  analysis  with  account  of 
apparent  masses.  For  a  spheroid  as  a  reference  body,  he  took  0.5  as 
a  baseline  values  for  both  coefficients  an  and  a33  and  varied  them 


within  the  shown  range.  Eighteen  years  later  Ludwig  and  Heins2 
increased  the  baseline  value  and  extended  the  range.  At  the  end  of 
the  1960s  based  on  experimentalresearch  by  Ibrahim,4'5  White  and 
Wolf 8  solved  the  five-DoF  equations  of  motion  using  approximately 
the  same  values  for  apparent  masses  and  introduced  apparent  mo¬ 
ment  of  inertia  a44.  Tory  and  Ayres12  implemented  theoretical  val¬ 
ues  computed  for  the  potential  flow  around  spheroid.  Not  only  did 
they  use  erroneous  equations,  but  their  modeling  data  seemed  to  not 
match  the  results  of  experiment. 

The  beginning  of  the  1980s  was  a  new  era  of  circular  parachute 
modeling.  Eaton13  and  other  research  teams14,16'17  developed  the 
right  set  of  six-DoF  equations  and  performed  stability  analysis. 
Eaton's  coefficients  if  referred  to  the  mass  of  a  he mispheroid rather 
than  a  spheroid  would  give  the  values  of  0.4  and  0.8  for  the  apparent 
masses,  which  is  fairly  close  to  those  of  Doherr  and  others.16'17 

In  spite  of  the  results  of  experimental  work  of  Yavuz  and 
Cockrell14  and  Yavuz15  showing  completely  different  scale  for  ap¬ 
parent  mass  terms,  it  looks  like  everybody  is  continuing  using  the 
classical  set  of  Doherr  (e.g.,  see  Refs.  18-20). 

The  present  study  was  not  aimed  at  tuning  these  coefficients, 
and  so  a  classical  set  of  apparent  mass  terms  was  used  here  also. 
However  the  authors  would  like  to  address  this  issue  in  the  future  .  It 
also  explains  why  element  a66  was  left  in  equations  notwithstanding 
that  for  undisturbed  canopy  it  is  always  neglected.  For  undisturbed 
canopy  it  contributes  to  the  damping  of  yaw  oscillations,  and  so 
authors  also  normalized  it  as  shown  in  Eqs.  (15)  for  the  future  study, 
letting  k66  however  to  be  zero  here. 

Computation  of  Forces  and  Moments 

Aerodynamic  Forces  and  Moments  for  Undisturbed  Canopy 

The  total  external  force  and  moment  acting  on  the  system  (Fig.  3) 
are  causedby  the  aerodynamic  effects  and  the  weight  of  each  system 
component.  Thus,  we  can  write 

F  =  Fa,d  +  G,  M  =  Majd  +  M‘a  (16) 

where  G  =  IRT [Q,  0,  mg]T  (the  apparent  (virtual)  masses  do  not 
contribute  to  the  weight  of  the  system). 

In  turn,  for  controlled  ADS 

Fa,d  =  F  +  F  M“ld  =M  +M-  (17) 

Let  us  address  further  in  this  subsection  aerodynamic  terms 
for  undistorted  canopy.  Initial  values  of  the  dimensionless  aerody¬ 
namic  coefficients  of  a  circular  parachute  were  obtained  from  two 
sources.  The  basic  shape  of  the  aerodynamic  curves  for  the  undis¬ 
torted  canopy  was  adopted  from  Knacke 26  whereas  initial  aerody¬ 
namic  coefficient  estimates  were  based  on  CFD  results  provided  by 
Mosseev.25 

The  functional  dependence  of  the  aerodynamic  coefficients  on 
the  angle  of  attack  obtained  by  Mosseev  is  shown  on  Fig.  6. 

On  this  figure  C D  (asp)  denotes  aerodynamic  drag  coefficient,  and 
C,„(as p)  denotes  total  aerodynamic  moment  coefficient,  both  de¬ 
pending  on  the  spatial  angle  of  attack. 

This  spatial  angle  of  attack  asp  and  its  components — angle  of 
attack  a  and  sideslip  angle  fi — are  shown  on  Fig.  7  and  can  be 
computed  using  {b  [-frame  components  of  an  airspeed  vector  V„  as 
follows: 


Table  2  Values  of  apparent  mass  coefficients 


No. 

Researcher 

Year 

Model 

Origin  of  [fe[ 

Ref.  body 

£ 

kn 

I'33 

A'44 

&15 

1 

Henn1 

1944 

Three  DoF 

O 

Spheroid 

0.5 

0.5  [0-0.6] 

0.5  [0-1.0] 

— 

2 

Ludwig,  Heins2 

1962 

Three  DoF 

O 

Spheroid 

0.5 

1.0  [0.6-1. 4] 

1.0  [0.6-1. 4] 

— 

3 

White,  Wolf8 

1968 

Five  DoF 

Sys.  c.g. 

Sphere 

1 

0.7 

0.7 

0.23 

4 

Tory,  Ayres12 

1977 

Six  DoF 

Can.  c.g. 

Spheroid 

0.5 

1.31 

2.12 

1.34 

5 

Eaton13 

1981 

Six  DoF 

Can.  c.g. 

Spheroid 

0.5 

0.2  [0-0.5] 

0.4  [0-1.0] 

Unclear 

6 

Yavuz,  Cockrell14'15 

1981 

Six  DoF 

Can.  c.g. 

Hemisphere 

1 

[1.31-6] 

[2.12-10] 

1. 5-3.0 

7 

Cockrell,  Doherr,  Saliaris 16,17 

1981 

Six  DoF 

O 

Hemisphere 

1 

0.5 

1.0 

0.24 

0.75 

8 

Present 

2003 

Six  DoF 

O 

Hemispheroid 

0.82 

0.5 

1.0 

0.24 

0.75 
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angle  of  attack  a  (°) 


Fig.  6  G- 12  parachute  aerodynamics. 


The  airspeed  vector  accounts  for  a  wind  vector  W„  in  the  follow¬ 
ing  manner: 

Va  =  V-ubRTWu  (20) 

Obviously  the  aerodynamic  force  vector  Fcanopy  depends  on  the 
spatial  angle  of  attack  and  dynamic  pressure  and  can  be  presented 
as  follows: 

Fcanopy  =  CD  (asp)qS0  (Va/\\Va\\)  (21) 

where  q  =  0.5p||  ||2  is  dynamic  pressure;  ||  •  ||  denotes  an  Euclidian 

norm  of  vector;  and  S0  =  ttRq  is  the  canopy’s  reference  area. 

White  and  Wolf8  showed  that  the  longitudinal  and  the  lateral  mo¬ 
tion  of  a  parachute  in  glide  plane  are  sufficiently  uncoupled.  There¬ 
fore,  the  longitudinal  and  lateral  motion  of  a  symmetric  parachute 
can  be  studied  separately,  similar  to  the  study  of  the  aircraft  lin¬ 
earized  dynamics.28,29  This  leads  us  to  assume  additionally  that  the 
longitudinal  and  the  lateral  motion  of  a  circular  parachute  in  a  glide 
plane  are  uncoupled. 

This  last  assumption  implies  that  roll  and  pitch  motion  of  the 
ADS  have  the  same  moment  characteristics,  that  is,  Croo  =  Cm  (yS), 
Cpitch  —  Cm  (a).  Therefore,  the  vector  of  aerodynamic  moment  can 
be  expressed  as  follows: 

Fro  11 

TF.'inopy  —  2qS0Ro  Fpuch  (22) 

_  C„  _ 


0 

▲ 

pma2 


Fig.  8  Yaw  rotation  scheme.  1  -4  PMA[  PMA:,^  0 


0  2  4  6  8  10 

time  (sec) 


Fig.  9  Effect  of  changing  riser  length. 

In  Eq.  (22)  C„  denotes  the  yaw  moment  caused  by  yaw  accelera¬ 
tion,  and  for  a  symmetric  body  it  is  equal  to  zero.  However  in  case 
of  transition  of  one  of  the  risers  from  one  state  to  another  while  one 
of  the  adjacent  risers  has  been  already  lengthened,  this  moment  is 
not  equal  to  zero.  The  next  subsection  addresses  this  issue. 

Yaw  Moment  Caused  by  Dynamic  Asymmetry 

As  just  mentioned,  significant  yaw  rotation  was  observed  during 
the  flight  tests  when  the  length  of  the  riser  was  changing  while  one 
of  the  adjacent  risers  has  been  already  lengthened.  Figure  8  clarifies 
this  situation. 

Here  “0”  stands  for  a  shortened  riser  ( lk  =  0),  and  “1”  denotes  a 
lengthened  riser  ( lk  =  1),  typically  while  riser  transition  (4 ...  5  s) 
yaw  angle  of  15  ...  20  deg  was  accrued  (Fig.  9). 

Therefore,  if  the  kth  riser  undergoes  transition  the  following  re¬ 
lation  is  valid: 

C„  =  sign(4_ !  -  Ik  +  i)C'  (/)/,  +  C'r  (23) 

In  Eq.  _  (23)  lk  =  (lk  -  -  (mm),  and  the  difference 

(l k  - 1  —  4  + 1 )  defines  the  sign  of  the  moment. 

By  analyzing  data  from  20  flight  tests,  the  functional  dependence 
C'  (/)  was  found  to  be  as  presented  on  Fig.  10  and  the  coefficient 
C'n  to  be  constant  and  equal  to  2  s. 

When  /  or  1  —  /  is  small,  as  happens  at  the  very  beginning  and  at 
the  very  end  of  the  actuation  process,  the  magnitude  of  the  rotation 
is  low  because  the  canopy  symmetry  distortion  is  minimal.  The 
rotation  is  at  its  highest  at  the  middle  of  actuation  process  when  the 
shape  distortion  is  maximal  (see  Fig.  9). 

Actuator  Forces  and  Moments 

Let  us  define  now  an  effect  of  riser  actuation. 

Vertigo’s  PMA27  used  as  a  prototype  of  control  riser  is  braided 
fiber  tubes  with  neoprene  inner  sleeves  that  can  be  pressurized.  In 
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Input 


Fig.  13  Block  diagram  of  the  system  identification  technique. 


Fig.  11  Dependence  /pma  =/fPpM  A  )• 


the  nominal  state  all  PMAs  are  pressurized.  Upon  venting,  PMA 
increases  its  length  by  approximately  one  third  from  /PMA  =  /min  to 
/pma  =  /max  =  IPma  +  A/pma  and  compresses  in  diameter.  [The  de- 
pendence/PMA  =  /(PPma)  where  -PPMA  stands  for  the  PMA  pressure 
is  shown  on  Fig.  11].  This  action  “deforms”  the  parachute  canopy 
creating  an  asymmetrical  shape,  essentially  shifting  the  center  of 
pressure  and  providing  a  drive  or  slip  condition  in  the  opposite  di¬ 
rection  of  the  control  action.  With  four  independently  controlled 
actuators,  two  of  which  can  be  activated  simultaneously,  eight  dis¬ 
tinct  control  actions  can  be  generated. 

The  change  in  the  aerodynamic  force  caused  by  PMA  activation 
^risers  was  modeled  as  a  function  of  PMA’s  relative  length  /,  number 
of  PMA  actuated  n,  and  involved  actuation  system  dynamics  with 
transition  time  r27: 

liters  (Oil  =f(J,n,T)  (24) 

Figure  12  shows  steady-state  values  of  this  dependence.  In  turn, 
the  actuator  moment  can  be  computed  as  Mrisers  —  OcP  x  Frisers, 
where  -PCp  =  LO,  0,  ~p]T . 

Model  Identification 

In  the  preceding  section  we  outlined  all  steps  involved  in  the  de¬ 
velopment  of  the  Six -DoF  model  for  a  controlledcircularparachute. 
Next  this  model  was  compared  with  the  flight-test  data.  This  com¬ 
parison  revealed  inadequacy  of  the  data  obtained  by  CFD  analysis 
(see  Fig.  6).  Therefore,  a  standard  nonlinear  system  identification 
technique32  was  applied  to  tune  the  CFD  dependences  CJ  (asp)  and 
C,*>sp)  (see  Fig.  13). 


First,  expressions  for  Cd  (<xsp)  and  C,„  (asp)  were  presented  in  the 
following  form: 

Cd(q:sp)  =  kCmCDa=0  +  kCoa  [CJ(asp)  —  CjjfO)]  (25) 

C,„  (<*sp)  =  kCm  C*  (asp)  (26) 

yielding  three  optimization  parameters  kCm,  kCDa ,  and  kCm.  The 
initial  value  of  Cd„=0  =  was  obtained  from  the  obvious 

equation  for  a  steady  descent  rate. 

Second,  an  appropriate  cost  function  was  chosen: 

r'f 

J=  /  IIP, (t)-Pu(t) ||  dt  (27) 

J  o 

where  Pu(t)  is  the  inertial  position  of  ADS  obtained  in  flight  test, 
P„(t)  denotes  the  estimated  position  of  ADS  obtainedin  simulation, 
and  tj  stands  for  the  flight  time. 

Therefore  the  objective  of  identification  was  to  minimize  cost 
function  (27)  by  varying  three  parameters  in  Eqs.  (25)  and  (26) 
(Fig.  13). 

Because  from  the  analysis  of  the  flight-test  data  it  became  obvious 
that  there  was  a  significant  difference  in  the  ADS  dynamics  during 
controlled  and  uncontrolled  drops,  the  identification  algorithm  was 
first  applied  to  the  data  obtained  from  the  uncontrolled  drop.  The 
resulting  values  of  kCm ,  kCDa ,  and  kCm  were  used  then  to  initialize  the 
second  step,  where  the  same  technique  was  applied  to  a  controlled 
drop.  Whereas  values  of  optimizationparametersobtainedin  the  first 
step  provided  estimates  of  ADS  aerodynamics  around  zero  angle 
of  attack,  their  adjustment  at  the  second  step  characterized  ADS 
dynamics  at  higher  angles  of  attack  with  nonzero  control  inputs.  As 
a  result,  some  averaged  values  were  found  to  match  both  controlled 
and  uncontrolled  real-drop  sets. 

All  two-  and  three-dimensional  trajectories  in  this  section  are 
plotted  in  East-North-Up  local  tangent  plane  (LTP)  coordinates. 

Uncontrolled  Model  Identification 

This  subsection  details  results  obtained  by  applying  the  identifi¬ 
cation  technique  to  the  uncontrolled  drop  data.  Because  the  model 
assumes  a  fully  deployed  canopy,  it  was  initialized  at  the  point  on 
the  real-drop  trajectory  where  the  canopy  was  fully  deployed. 

Figures  14-16  contain  the  three-  and  two-dimensional  plots  of 
the  trajectories,  obtained  in  one  of  the  real  flight  three-DoFand  six- 
DoF  simulations.  (The  three-DoF  model  of  the  parachute  was  used 
in  the  preliminary  study.21) 

Note  that  the  real  trajectory  has  three  points  where  the  parachute ’s 
direction  of  motion  changes  more  than  90  deg.  This  particular  data 
set  was  selected  because  of  the  richness  of  the  frequency  spectrum 
of  the  parachute  motion  involved.  Fig.  17  contains  the  wind  profile 
for  the  same  drop.  It  clearly  indicates  that  both  wind  direction  and 
its  magnitude  experience  significant  changes  with  time. 

Because  during  the  uncontrolled  drop  the  angle  of  attack  and 
sideslip  angle  are  close  to  zero  the  value  of  kCm  was  set  to  zero. 
This  choice  implies  that  the  amplitude  of  coning  motion  can  differ 
from  the  flight-test  data;  however,  its  natural  frequency  does  not 
depend  on  kCm  for  uncontrolled  drops.  Finally,  the  optimal  values 
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Fig.  15  Horizontal  projection. 


kCd0  and  kCDa  were  obtainedby  determininga  priori  the  range  of  their 
possible  values  from  physical  considerations  and  applying  simple 
search  techniques  over  the  resulting  grid  of  this  parameter  space. 

Figures  18-20  include  plots  of  the  x,  y,  and  z  components  of  the 
inertial  velocity  obtained  from  the  six-DoF  model  simulation  and 
during  the  flight  test.  Analogous  data  are  plotted  for  the  three-DoF 
model.  A  conspicuous  feature  of  the  six-DoF  model  can  be  seen  in 
the  plots  of  the  vertical  or  z  component.  (Figure  20  shows  absolute 
values).  Although  the  vertical  component  of  the  three-DoF  model 
remains  constant  during  the  flight,  the  same  component  of  the  six- 
DoF  model  follows  the  flight-test  data  reasonably  well.  Therefore 
optimization  over  a  set  of  two  parameters  gave  fairly  good  results. 
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Fig.  17  Wind  profile. 


Figure  21  contains  power  spectrum  density  (PSD)  plots  for  Euler 
angles  (roll,  pitch,  and  yaw)  taken  from  the  flight  test  and  the  six- 
DoF  model  run.  Among  dozens  of  methods  available,  the  Multiple 
Signal  Classification  method  developed  by  Schmidt33  was  applied. 
This  method  was  chosen  to  demonstratea  general  similarity  between 
real  and  simulated  data  as  well  as  to  emphasize  the  matching  of 
principal  eigenfrequencieson  the  both  sets  of  data. 

Numerical  values  of  these  principal  eigenfrequencies  for  each 
channel  are  depicted  in  Table  3. 

Similar  results  for  the  principal  eigenfrequencies  were  obtained 
independently  for  another  set  of  ADS  flight-test  data  using  Welch’s 
method  as  well.34  That  shows  the  consistency  of  ADS  flight-test 
data  taken  in  the  presence  of  the  different  wind  profiles,  indicating 
that  the  data  in  Table  3  reflect  the  inherent  motion  of  the  ADS,  and 
not  the  wind  spectrum. 

In  turn,  the  wind  spectrum  characteristics  for  the  real  drop  being 
analyzed  are  presented  in  Fig.  22.  The  results  clearly  show  that 
the  wind  energy  was  concentrated  around  much  higher  frequency 
(1.9  Hz)  rather  than  0.1 1 -Hz  parachute’s  harmonic  just  seen. 

The  main  conclusion  to  be  drawn  from  spectral  analysis  is  that 
the  six-DoF  model  dynamics  are  sufficiently  close  to  those  of  the 
actual  system. 
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Controlled-Drop  Dynamics  Verification 

To  isolate  the  influence  of  control  on  the  parachute  aerodynamics, 
the  following  approach  was  used.  First  the  wind  profile  measured 
during  the  controlled  drop  was  used  to  simulate  the  uncontrolled 
model.  Then  this  simulated  data  were  subtracted  from  the  flight- 
test  data,  and  the  result  was  used  as  an  input  data  for  the  system 
identification  algorithm.  This  looks  reasonable  because  the  con¬ 
trol's  influence  on  the  descent  rate  is  minor;  hence,  the  new  input 
data  are  supposed  to  approximate  controlled  drop  in  the  absence 
of  wind. 

An  example  of  trajectories  at  the  presence  of  strong  wind  is 
shown  in  Fig.  23.  By  applying  the  described  procedure,  a  con¬ 
tribution  of  control  actions  was  received  as  a  difference  of  two 
trajectories. 

The  result  of  the  identification  procedure  for  this  case  is  depicted 
in  Fig.  24.  Here  the  “Control  +  NoWind”  graph  represents  the  posi¬ 
tion  obtained  by  driving  the  six-DoF  model  with  no  wind  and  with 
control  inputs  recorded  in  the  real  drop.  (The  logic  of  the  control 
algorithm  is  described  in  detail  in  Ref.  35.) 

Finally,  Fig.  25  compares  controlled-drop  trajectories  obtained 
in  simulation  with  double  tuned  aerodynamics  and  flight  test. 


Table  3  Eigenfrequencies 


Euler 

angle 

PSD  rate,  Hz 

Real  drop 

Six-DoF  model 

Roll 

0.109 

0.114 

Pitch 

0.112 

0.111 

Yaw 

0.12,  0.22,  and  0.406 

0.13,  0.22,  and  0.41 

Table  4  Values  of  optimization  parameters 


Optimization 

parameter 

Initial 

value 

Uncontrolled 

drop 

Controlled 

drop 

kcDO 

i 

1.2 

1.2 

kCDa 

i 

1.005 

1.25 

kcm 

i 

0 

0.6 

0  0.1  0.2  0.3  0.4  0.5  0.6 

frequency  (Hi) 

Fig.  21  Euler  angles  PSD  analysis. 


Fig.  20  Comparison  of  absolute  values  of  Vz  components. 


Fig.  23  On  the  control  contribution. 
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Fig.  24  Control  contribution  comparison. 


Clearly  results  obtained  in  simulation  approximate  flight-test  data 
very  well. 

Final  values  of  optimization  parameters  are  shown  in  Table  4. 

Conclusions 

The  paperpresentsresults  of  the  development  of  a  six-DoF  model 
of  the  controlled  circular  parachute.  Key  contributions  of  the  paper 
include  the  following:  1)  controlled  model  of  the  six-DoF  circular 
parachute  based  on  detailed  modeling  of  the  ADS  geometry  and 
utilization  of  the  actuators  performance  data  provided  by  the  man¬ 
ufacturer;  and  2)  successful  application  of  the  two-step  nonlinear 
system  identification  algorithm  to  refine  the  analytical  values  of  the 
aerodynamic  coefficients  based  on  available  flight-test  data. 

The  resulting  model  acquitted  itself  well  when  compared  with 
the  flight-test  data.  However,  further  improvement  is  possible.  It 
includes  involving  more  optimization  parameters  into  the  model 
identification  procedure.  For  instance,  rather  than  finding  averaged 
values  of  the  aerodynamic  force  coefficients  the  effect  of  zero,  one, 
and  two  risers  activation  might  be  addressed  separately.  Optimiza¬ 
tion  parameters  can  also  cover  coefficients  of  the  apparent  mass 
tensor  terms  as  a  function  of  the  angle  of_attack  and  canopy  shape 
while  applying  control  inputs  kmn  =  kmn  (/,  n).  More  complex  cost 
function  (adequacy  criterion)  or  even  multiple  criteria  might  be  con¬ 
sidered  also.  Further  analysis  could  also  include  controlled  ADS 
asymmetry,  meaning  the  accounting  of  more  terms  in  the  inertia 
and  apparent  mass  tensors. 

Appendix:  Geometry  and  Masses  of  the  G-12-Based 
Aerial  Delivery  System 

Tables  A1-A3  contain  geometric  and  mass  data  used  to  estimate 
the  position  of  static  mass  center,  the  moments  of  inertia  (Table  A4), 
and  arms  of  corresponding  forces  while  computing  moments  acting 
on  the  generic  G-  12-based  ADS. 


Table  A1  ADS  components’  dimensions 


Component  Value 


Length  of  suspension  lines  /sl,  m  15.55 

Length  of  pressurized  PMA  /pma.  m  5.80 

Length  of  vented  PMA /*MA,  m  7.62 

Radius  of  uninflated  canopy  Rq,  m  9.75 

Radius  of  inflated  canopy  Rp  =  jRq,  m  6.50 

Canopy  shape  ratio  s  0.82 

Cone  half-angle  y ,  deg  15.31 


Dimension  of  a  cubic  payload  container  ap\,  m  1 .22 


Table  A2  ADS  components’  masses 


i 

Component 

Mass,  kg 

i 

Canopy  m i 

22.8 

2 

Suspension  lines  m2 

35.3 

3 

Four  PMAs  m3 

13.4 

4 

Payload  m 4 

990.0 

ADS  total 

1061.5 

Table  A3  Auxiliary  geometric  relations 


z  coordinate  of  a  specific  point 

Value,  m 

Centroid  of  hemispheroidal  canopy  z\ 

-2.76 

Centroid  of  suspension  lines  Z2 

7.50 

Centroid  of  PMAs  c3 

17.80 

Centroid  of  payload  Z4 

21.20 

Static  mass  center  of  the  whole  ADS  zq 

20.19 

Center  of  pressure  zp  =  —^sRp 

-2.00 

Table  A4  Moments  of  inertia,  kg  •  m2 


i 

Component 

^aa  ~  ^bb 

4  =  4 

4  =  4 

1 

Canopy 

365.93 

539.16 

623.98 

2 

SLs 

1,047.09 

3,032.22 

770.76 

3 

PMAs 

53.43 

4,296.81 

36.98 

4 

Payload 

245.59 

445,287.93 

245.59 

ADS  total 

— 

453,156.13 

1,677.30 
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